Detecting maximally entangled states without making the Schmidt decomposition 
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The bipartite entanglement of a pure quantum state is known to be characterized by its Schmidt 
decomposition. In particular the state is maximally entangled when all the Schmidt coefficients are 
equal. We point out a convenient method which always yields a single analytical condition for the 
state to be maximally entangled, in terms of its expansion coefficients in any basis. The method 
works even when the Schmidt coefficients cannot be calculated analytically, and does not require 
their calculation. As an example this technique is used to derive the Bell basis for a system of two 
qubits. In a second example the technique shows a particular state to never be maximally entangled, 
a general conclusion that cannot be reached using the Schmidt decomposition. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.67.Hk, 02.10.Ud 



I. INTRODUCTION 

Entanglement or the existence of quantum correlations 
between physical systems is currently of great interest 
both theoretically and experimentally, see Ref. [l| and 
references therein. From the theoretical point of view 
the study of entanglement is shedding light on the ba- 
sic nature of quantum mechanics Q , and providing sur- 
prising connections to other areas, such as the study of 
black holes 0]. From a practical point of view it has 
been realized that entanglement is a resource that can 
enable information processing tasks that turn out to be 
impossible or inefficient when tackled by classical ma- 
chines [3]. These capabilities include tclcportation, cryp- 
tography, and computation [l| . 

In this article we will consider the entanglement of 
pure bipartite states, i.e. states which describe globally 
pure quantum mechanical systems which have been par- 
titioned into two subsystems. As is known, the Schmidt 
decomposition of a pure bipartite stateprovides a qual- 
itative measure of its entanglement P, Q ■ On the other 
hand the von Neumann entropy of either subsystem pro- 
vides a unique quantitative measure of the degree of en- 
tanglement of the whole quantum state [l|, Q • 

In particular we will consider maximally entangled 
states i.e. those states of a quantum system which con- 
tain the largest amount of entanglement possible. Since 
entanglement is a resource for information processing, 
these states are of particular importance. Indeed pro- 
tocols for quantum key distribution, dense coding and 
teleportation rely on maximally entangled states In 
view of the inevitable noise and decohcrence that accom- 
panies an experiment, maximally entangled states also 
provide a benchmark for other, less entangled, states, 
as in the process of distillation The Bell states are 
perhaps the best known example of maximally entangled 
bipartite pure states. 

How do we know if a pure bipartite state is maximally 
entangled ? We can decide using either of the measures 
mentioned above. If the state is entangled maximally all 
the coefficients in the Schmidt decomposition are equal; 



also the reduced von Neumann entropy is maximised to 
its upper bound log d, where d is the dimension of the 
smaller partition. However if d > 4 neither of these mea- 
sures can provide analytical conditions (on the expansion 
coefficients of the state in some basis), and have to be 
computed numerically. 

In this article we point out an existing technique from 
algebra that can provide a single analytical condition (for 
arbitrary but finite d) that the coefficients (in any basis) 
of a maximally entangled state have to obey. The method 
therefore is capable of deciding quite generally if a certain 
state is maximally entangled or not. We demonstrate 
below how this method shows a certain state can never 
be maximally entangled, a general conclusion that cannot 
be reached using the Schmidt decomposition or the von 
Neumann entropy. In addition the method can also be 
used to obtain the maximally entangled basis given the 
dimensions of the two subsystems. As an example below 
we use the method to derive the Bell basis for a system 
of two qubits. 

The rest of the article is arranged as follows. In Sec- 
tion[n]we outline the general algorithm. In Section HTTl we 
use the technique to derive the Bell basis for a system of 
two qubits providing also a detailed commentary on our 
method, and its simple implementation in Mathematica. 
In Section IIVI we consider a case which cannot be solved 
analytically using the Schmidt decomposition or the von 
Neumannn entropy. Section IVl provides a Conclusion. 



II. GENERAL FORMALISM 

The steps of our technique are as follows: 

1. Begin with a quantum state lip). 

2. Find the density matrix p = 

3. Trace over the states of one subsystem, say A, to 
find the reduced density matrix ps ~ Tr^p. 

4. Find the characteristic polynomial P[pB,x] of this 
matrix with respect to a dummy variable x. 
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5. Find the subdiscriminant sequence Dq[P], of the 
polynomial P, where q = 1, 2 . . . d Q. 

6. If the last but one member of the sequence, i.e. 
Dd~i[P], equals zero then lip) is maximally entan- 
gled. 

The basic tool we use is the subdiscriminant sequence 
of the characteristic polynomial of the reduced density 
matrix of cither subsystem. The subdiscriminant se- 
quence of any polynomial can be found in textbooks on 
algebraic geometry such as Ref. @. A detailed intro- 
duction for physicists has been provided in Ref. and 
will not be repeated here, although we will discuss a few 
important points. In general the sequence contains d 
members. The various members denote the number of 
repeating zeros of the polynomial. For example the first 
member of the sequence is the discriminant of the poly- 
nomial P 



Di[P]=l[{X., 



(1) 



where Ai are the roots of the polynomial P. Since density 
matrices are Hermitian, the Xi are all real. From Eq. ^ 
we can see that vanishes whenever two or more 

roots of P are equal. 

For the purpose of this article the most important 
member of the subdiscriminant sequence is its second to 
last entry 



i?,_i[P]=^(A,-A,)2 



(2) 



which equals zero only if all the eigenvalues of the re- 
duced density matrix are equal. In turn this implies 
the equality of the Schmidt coefficients and maximal en- 
tanglement of the state. It is an important fact that 
£'i[P], Dd-i[P], and actually the whole subdiscriminant 
sequence Dq [P] can always be found analytically in terms 
of the coefficients of the polynomial P [g| . Significantly 
this does not require calculation of the eigenvalues Xi, i.e. 
one need not make the Schmidt decomposition 

Generally i5d_i[P], which has to vanish for a maxi- 
mally entangled state, is itself a polynomial in the co- 
efficients of expansion of the given quantum state writ- 
ten in an arbitrary basis. It follows that from a broader 
perspective the present technique maps the problem of 
finding maximally entangled states of a bipartite system 
to that of finding the roots of a multivariate polynomial. 
In order to expose the working details of our procedure 
we present two examples. We first rederive the familiar 
Bell states using our method; then we consider a more in- 
volved example which cannot be solved analytically using 
Schmidt decomposition or the von Neumann entropy. 



III. A 2 X 2 SYSTEM : THE BELL BASIS 

We consider a system which is divided into two subsys- 
tems A and B. Each subsystem contains a qubit which 
can exist in a superposition of the states |0) and |1). An 
arbitrary unnormalized pure quantum state of this sys- 
tem can be written as 



1^) ==p|00) +g|ll) +r|10) +s|01). 



(3) 



Here we have used the product basis, which is often a 
convenient one to use, although for the method to be 
demonstrated \^) can be expressed in any basis. For the 
particular task of deriving the Bell states we have chosen 
the coefficients (p, g, r, s) to be all real. 

The density matrix corresponding to the state in 
Eq. ([2]) can be written easily and its trace over the first 
qubit yields the reduced density matrix 



^ ' ps + qr q^ + 
where we have used the representation 



(4) 



1 





1 



,|0)(1| 
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(5) 

etc. The characteristic polynomial oi pb [Eq- O] is 

P[pB,x\ =x^ - {p^ + q^ +r'^ + s^)x+ {pq - rs)^ (6) 

where a; is a dummy variable. Clearly P is quadratic 
in X and has two roots which are non-zero in general. 
They are the eigenvalues oi pb and yield the Schmidt 
coefficients 0; as promised, we will not calculate them. 

The Subresultants[P, P', x] function in Mathemat- 
ica directly yields the subdiscriminant sequence with the 
polynomial P and its derivative P' with respect to x as 
inputs. Wc note that Mathematica requires that the co- 
efficient of the highest power of x in P as well as P' to be 
1. This can be arranged easily. Also in the general case 
the number of terms in the sequence equals the degree of 
P, i.e. d. In the present example therefore there are two 
terms in the sequence. The penultimate term in the sub- 
discriminant sequence of P in Eq. © is the discriminant 
of P, Di[P], which equals zero whenever the two roots 
of Pb coincide Q . We find 

D,[P] = [{p + qf + (r - sf] [{p - qf + (r + sf] . 

(7) 

If It/i) is to be maximally entangled -Di[P] = 0. We can 
extract the Bell basis using this criterion. Specifically, 
we can see that {p ^ q = Q,r = ±s) is & solution set and 
yields, using Eq. ([3]), the states 



1^ 



1,V 



|io)±|oi) 

^/2 ■ 



(8) 



after normalization. Similarly, (p — ±q, r = s = 0, ) is a 
solution set and yields the states 



1^-3.4) 



|00) ± 111) 
V2 ^ 



(9) 
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after normalization. As is well known, the states \'4'i.2,3a) 
constitute the Bell basis Any state, and therefore any 
maximally entangled state can be expressed in this basis. 
For example, the solution set [p = —q,r = s) yields the 
(unnormalized) maximally entangled state 



V2\iJ4)+rV2\iJi). 



(10) 



We note that (p, g, r, s) can be complex in general; how- 
ever the analysis in that case is not very different from 
that presented here (also see below). 



IV. A 5 X 5 SYSTEM 

The superiority of the method proposed in this article 
over that of Schmidt decomposition or von Neumann en- 
tropy calculation becomes clearer for systems of higher 
dimension, i.e. qudits. Physically qudits are of interest as 
some of the proposed candidates for quantum computa- 
tion possess d > 2, such as a molecule with ro-vibrational 
states 01 and an alkali atom with hyperfine-Zecman lev- 
els [|. 

Here we consider a system with two parts each con- 
taining a 5-level system with states |0), |1), |2), |3), and 
1 4). Specifically, we consider the unnormalized state 

IV') = |02) + 2|10) + |20) + |21> + |22> + |23) 

+ |24)+3|33)+p|44), (11) 

where p is the only unknown real coefficient of expansion 
in the product basis. Tracing the density matrix over the 
five levels of system A, we obtain 



/ 5 1 1 
'ill 



PB 



1 1 

1 1 

112 1 1 

1 1 1 10 1 



Villi 



p- 



1/ 



(12) 



where we have used a matrix representation analogous to 
Eq. ([5]). The characteristic polynomial of [Eq. p2)) ] 
is 

P[pB,x] = -x^ -f + 19) - (I8p2 4- 105) + 

(91p2 + 183) x^ - (134/ + 72)x + 36/.(13) 

P is a quintic, and not generally solvable analytically 
in terms of radicals [H]- This implies that the Schmidt 
coefficients have to be found numerically. However the 
subdiscriminant sequence can be found analytically and 
its fourth member, apart from an irrelevant numerical 
prefactor, is 



is therefore maximally entangled when D4[P] 
which yields the four roots 



1/2 



1/2 



p = ± 



, ± 



(15) 



2 / ' \ 2 / 

Since none of these solutions are real, \ip) can never be 
maximally entangled. 

We note that if we had assumed p to be complex in 
Eq. (ini), then Eqs. (HH), and HH) would undergo 

the transform p^ — > The conditions presented in 

Eq. (fTS]) would then hold for |p| and would be impossi- 
ble to achieve since by definition the modulus of p, 
is real. Therefore even if p is allowed to be complex lip) 
can never be entangled maximally. We also note that the 
case of Eq. (fT4|) is exceptional in that it is an analytically 
solvable equation, i.e. a quartic, in the variable of inter- 
est, p. This is atypical. Although the subdiscriminant 
sequence can always be obtained analytically, the roots 
of its members typically have to be found numerically. 

The above demonstration, although it uses a somewhat 
arbitrary quantum state, shows the general usefulness of 
the method introduced in this article. The inability to 
maximally entangle the state of Eq. pT|) for any value 
of p cannot be established generally by calculating the 
Schmidt coefficients or the von Neumann entropy. 



CONCLUSION 



In this article we have pointed out a technique that 
always yields a single analytical condition that is satisfied 
by the coefficients of a pure bipartite quantum state if it 
is maximally entangled. The method is superior to that 
of Schmidt decomposition and von Neumann entropy for 
multilevel systems of dimension greater than 4. Further, 
it can also be used to obtain the maximally entangled 
basis for systems of arbitrary but finite dimension. 



D4[P] = 2p^ - 14p2 + 197. 



(14) 
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